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Abstract 

We study the growth of a periodic pattern in one dimension for a model of 
spinodal decomposition, the Cahn-Hilliard equation. We particularly focus on 
the intermediate region, where the non-linearity cannot be negected anymore, 
and before the coalescence dominates. The dynamics is captured through the 
standard technique of a solubility condition performed over a particular fam- 
ily of quasi-static solutions. The main result is that the dynamics along this 
particular class of solutions can be expressed in terms of a simple ordinary 
differential equation. The density profile of the stationary regime found at 
the end of the non-linear growth is also well characterized. Numerical simula- 
tions correspond satisfactorily to the analytical results through three different 
methods and asymptotic dynamics are well recovered, even far from the region 
where the approximations hold. 
Pacs numbers : 05.45.Yv, 47.20. Ky, 47.54. +r 
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I. INTRODUCTION 



When a homogenous system departs suddenly from equihbrium, the fluctuations around 
the initial ground state are linearly amplifled and for example the homogenous phase sepa- 
rates spontaneously into two different more stable states. The interfaces which delimit the 
domains of each phase form a complex pattern and interact with each other, giving rise to 
interface dynamics or pattern formation. Its results can be a slow process of coarsening 
which ends up with only two well separated domains. This process of flrst order phase 
transition arises particularly for binary mixtures or alloys vapor condensation [Q, 
ferromagnetic Ising model or thin fllms of copolymers 0] . 

For the most general, flrst order transitions initiate in two different ways: flrst, a nucle- 
ation process, where the homogenous state is put suddenly in a metastable conflguration, 
and an energy barrier has to be crossed before the transition appears. This is the typical 
dynamics of cavitation, for instance 0. The other method is spinodal decomposition where 
the system leads in a linearly unstable conflguration; such is the situation that we will study 
here. In this latter case, three different regimes are identifled in the dynamics: flrst the 
linear instability of the homogenous phase develops from the fluctuations, leading to the 
creation of a modulation of the order parameter at a well deflned length scale. The mod- 
ulations grow exponentially with time as long as the non-linearities are negligeable. This 
stage is very short and results mainly in the selection of a particular length scale for the 
process. Non-linearities rapidly slow down the growth of the modulation resulting in an 
interface pattern composed of well deflned interfaces delimiting domains containing one of 
the two stable phases. Remarkably, this intermediate stage conserves quite perfectly the 
modulation width, so that the resulting pattern is of almost the same length scale than the 
one selected initially. Finally, a slow, self-inhibiting dynamics dominates the last stage of 
the process, due to the interactions between the interfaces. The different regions of each 
phase coalesce in the so-called Ostwald ripening where the number of domains diminishes 
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whereas their typical size increases. The asymptotic state is decomposed into two domains, 
one for each phase. This coarsening dynamics is in fact present aheady from the beginning 
of the spinodal decomposition; however, as we will discuss below, its influence on the two 
first stages can be often neglected. 



Hillert Cahn and Hilliard have proposed a model equation for a scalar order 
parameter describing the segregation for a binary mixture. This model, known as the Cahn- 
Hilliard equation (C-H later on), belongs to the Model B class in Hohenberg and Halperin's 
classification Indeed, different models of phase separation have been proposed, depend- 
ing whether the order parameter is a scalar or a vector, or wether it is or is not a conserved 
quantity (for a review see The (C-H) equation is in fact a standard model for 

phase transition with conserved quantities and has applications to phase transition in liquid 
crystals [|10|, segregation of granular mixtures in a rotating drum [^, or formation of sand 
ripples |r2| , p!3| . It is a partial differential equation to which a conservative noise is added to 



account for thermal fluctuations 



Figure (|I|) shows snapshots of the numerical simulation of the (C-H) dynamics which 
represents the full phase transition process after a quench in temperature. In that case, 
thermal fluctuations have been omitted in the dynamics, but were present in the initial 
conditions. The three main stages of the spinodal decomposition described above are clearly 
distinguished: first, from Fig. (|I] (a)) to Fig. ([l| (b)), we observe the selection of a typical 
length scale for the modulations, then the non-linear growth and its saturation from Fig. (|l| 
(b)) to Fig. (|l| (c)). We note that the number of peaks has been almost conserved between 
these two configuration; on the other hand the amplitude of the modulation has now reached 
almost its asymptotic value and will not change significantly in the further dynamics. On 
the contrary, the coarsening dynamics is observed between Fig. (c)) and Fig. ([1| (d)) and 
the typical length of the pattern is increasing. 
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FIGURES 




FIG. 1. Time evolution of the order parameter ^{x,t) for e = —1, dx = 0.1227. (a) initial 
conditions at t = are taken randomly with a very low amplitude (5 • 10~^) ; (b) at time t = 15, 
the amplitude of the modulation has decreased, while only long wavelength contributions are still 
present. The small scale perturbations have been damped by the (C-H) dynamics ; (c) at t = 225, 
the modulation has almost reached its final amplitude, keeping roughly the same number of peaks 
as before ; (d) at t = 1800, we observe that the number of domains has decreased from the 
coarsening dynamics. 



Note that in the numerical experiments, contrary to a realistic experience, we go instan- 
taneously from one temperature to an other (ideal quench). And we will see in the following 
that the system will evolve in a smoother ways in the sense that, due to finite diffusion 
process, it will not instantaneously reach the thermodynamical state associated with the 
temperature of the quench. 

In this paper, we will focus on the one dimensional (C-H) equation and our aim is to 
offer a consistent description for the pattern formation, corresponding to the intermediate 
(non-linear) regime. We will obtain a simple ordinary differential equation describing the 
dynamics along a family of quasi-static periodic solutions. We recover the linear regime for 
short times, and correctly reproduce the saturation of the second (non-linear) stage, in the 
case of small initial perturbations, in the region close to the critical point (i.e., for a symmet- 
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ric mixture). These results are valid in the limit where thermal fluctuations and coarsening 
processes are both neglected. We will discuss these important assumptions and show how 
the resulting ordinary equation depends on the wavelength of the periodic solution. 



As mentionned above, numerous models for phase transition have been proposed ; an 
important activity has been devoted to the description of their dynamics, using both statis- 
tical methods or numerical simulations (for a review see [|15|). However, these works mainly 
concentrate on the late stage of the spinodal decomposition where the coarsening dynam- 
ics dominates and which exhibits " dynamical scaling" : the dynamics presents a self-similar 
evolution where time enters only through a length scale L{t), associated with a typical length 
of the domains or the rate of decay of the inhomogeneities. For instance, scaling arguments 
and stability criteria give the law L(t) ~ t^^'^ for spatial dimensions greater than one and a 
logarithmic behavior for one dimension in the case of the (C-H) equation . 



In this article, by spinodal decomposition, we refer to the two first stage only, excluding 
the coarsening dynamics or the third stage. 

Only little is known experimentally of these two regimes of the dynamics: indeed, they 
are too brief and therefore very hard to capture. However, these stages were observed in a 
recent experiment on two-dimensional demixion of copolymer @] which motivated our work 
since it shows the need for a better understanding of the dynamics before the coalescence. 
While the linearized theory gives a full understanding of the first stage, the second stage 
of phase separation, which concerns the saturation of the growth through the non linearity, 
appears to have been less studied. There exists numerical attempts to provide descriptions 
of the saturation of the profile up to its stationary regime, using for example a concentration 
dependant diffusion coefficient, ||T6[ D = — $), which leads to a modified Kuramoto- 

Shivashinski equation and enables to have a saturation of growth. Here, on the contrary, we 
work with a constant diffusion coefficient : the non-linearity will only come from the usual 

term of the Landau free energy. 

The paper is organized as followed: first, we present a brief review of general properties 



of phase segregations and on the (C-H) model, mainly to fix the notation. We will reproduce 
briefly the original derivation by Cahn and Hilliard, and we will restrict ourselves to the one 
dimensional case. In part III, we motivate the different assumptions of our calculations. 
Numerical simulations are used to determine the role of the noise and the influence of the 
coarsening in the early dynamics. Then, in part IV we focus on interfaces ; in particular, we 
will exhibit a two parameters family of solutions, speciflc of the one dimensional case, the 
so-called soliton lattice. Finally in section V, we will make use of the solvability criterion 
in order to select the dynamical evolution of the density proflle among a selected "ansatz" 
solution. Eventually, we compare these results with numerical studies of the full (C-H) 
dynamics shown at the end. We conclude with a discussion of possible extensions of this 
work. 

II. THE CAHN-HILLIARD MODEL 

The Cahn-Hilliard theory is a modifled diffusion equation ; it is a continuous model, 
which reads in its dimensionless form: 



Here r and t represent the position vector and the time, the vectors being noted with 
bold fonts. $ is the order parameter, a real number ; for instance, it can correspond to 
the dimensionless magnetization in Ising ferro magnet, to the fluctuation of density of a 
fluid around its mean value during a phase separation or to the concentration of one of 
the components of a binary solution in some region around Y^. D is the diffusion constant 
and e is the dimensionless control parameter of the system ; it is often identifled to the 
reduced temperature {e = '^-^ where Tc is the critical temperature of the phase transition). 
This equation, first derived by Cahn and Hilliard [0, has also been retrieved by Langer 
from microscopic considerations. As written, the (C-H) equation does not account 
for thermal fiuctuations present in the system. They can be added through a Langevin 




(1) 
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force, which integrates in the Fokker-Planck equation for the probabihty distribution of 
||17|| . However, as explained by the thermal fluctuations can equivalently be taken 



in account through a random noise term on the r.h.s. of equation (|I|). Thus, the (C-H) 
equation reads, in its more general form: 

^ (r, t) = V2(|$ + 2$^ - + C(r, t)), (2) 

where C is a white noise of norm unity, whose amplitude is proportional to the square 
root of the temperature of the system. 

The (C-H) model is a conservative model for the order parameter $. Indeed, it can be 
written as: 

where j is the current associated with $. Moreover, this current obeys the standard law 
related to the gradient of a so-called chemical potential (j = — V/i). For (C-H), /i is itself 
defined as the functional derivative of a free energy F, through: 



with F being in that case the usual Landau- Ginzburg density: 



F=^ (^(V$)2 + + $4 

The homogeneous stationary solutions for the noiseless (C-H) equation are extrema of the 
effective potential V{^) = e^"^ + <l>'^. For positive e, there is only one homogenous solution 
$ = which is linearly stable ; for negative e, the stationary solution $ = undergoes a 
pitchfork bifurcation and three stationary solutions exist. $ = is still a stationary solution, 
but it is now linearly unstable ; two other symmetric solutions $ = ±-^^ are stable and 
have the same free energy F = — e^/32. 



Thus, a first order transition can be experienced by quenching the system suddenly from 
a positive reduced temperature e to a negative one. Spinodal decomposition is the resulting 
dynamics. Since for all positive e the system is described by $ = 0, we only have to study 
the case where we start at t = with $ = and a negative e. This is what was shown on 
figure (|I]) in one space dimension, where the noise has been omitted except for the initial 
condition, where it consists of a random noise of small amplitude around the mean value 
< $ >= 0. This can be justified since the noise level, being proportional to the square 
root of the temperature, is higher before the quench than after. Thus, taking a noisy initial 
condition and omitting the noise further on can be interpreted as neglecting the noise of the 
quenched system compared to the residual noise coming from the "hot" initial temperature. 
However we will discuss more precisely below the influence of the noise in the quenched 
phase. 

When the equation is studied for a constant £, via a rescaling of $ (as v^— position 
r (as r/ yj—e) and time (as we observe that we could restrict the dynamics to the 

case e = —1. However, since we will later on compare stationary solutions of (C-H) with 
different reduced temperature, we will continue to write the equation with a given keeping 
in mind that the dynamics can always be rescaled to the case e = —1. 

The stability of the solution $ = can be studied by linearizing equation ([l|) around 
$ = (i.e. neglecting the non linear term $^); considering $ as a sum of Fourier modes: 

$(r,t) = 5:0qe^''-+'^* 
q 

where (pq is the Fourier coefficient at t = 0, we obtain that the amplification factor a{q) 
satisfies: 

a(q) = -(g2 + £)g2 
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It shows immediately that $ = is hnearly stable for e > while a band of Fourier 
modes are unstable for negative e, since (j(q) > for < g < {—e/2). Moreover, the most 
unstable mode (where a is maximal) is for = ^— £:/2(with am = fg)- We can anticipate 
that this wave number of maximum amplification factor will dominate the first stage of the 
dynamics; in particular, it explains why the modulations appear at length scales close to 
Xm = 27T/qm, tlic wavc length associated with g^- Indeed, we show on figure (H) the time 
evolution of the usual structure factor in one dimension: 

5(g) = <l'(g)<l>(g)* 

where $ is the Fourier transform of the field $ ($* standing for its complex conjugate). 
We have taken the noiseless (C-H) equation with random initial conditions; the curve is 
obtained through an average over 100 initial conditions. 




0.01 0.1 

wave number k 
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FIG. 2. The structure factor S{k) as a function of the wave number k for different time 
t = 0, 50, 100, 150 and t = 500 time units; the higher the peaks are, the larger the time. The 
curves are an average over 100 initial conditions taken as a random noise of amplitude 5 • 10~^, 
the discretization is over 4096 grid points with the grid space dx = 0.6, and e = —1. At t = 0, we 
observe the flat spectrum of the white noise. For t = 0.5, the spectrum reflects the amplification 
fator: the peaks of the factor is located at qm = 0.5 while for all the modes q > \f2 the initial noise 
has been damped. Then for t = 100 and t = 150 we observe the formation of higher harmonics 
but the peaks of the structure factor stay around qm- However, at t = 500. the coarsening of the 
solution has began since the maximum of S{k) is now at a larger wave length. 

The different regimes are again well identified: at short times we see that the modulations 
whose wave number is close to qm grow rapidly from the white noise, while the fluctuations 
for q > V2qm for which the amplification factor is negative are damped. Then, higher wave 
numbers emerge, which correspond roughly to harmonics mode of the initial modulations. 
It corresponds to the intermediate stage of the dynamics, where the single mode approxima- 
tion of the profiles is not valid any more and the dynamics is in a highly non-linear regime. 
Notice however that the structure factor keeps its peak located around qm] as we will dis- 
cuss below, it is indicating that the number of domains stays almost unchanged during this 
regime. Later on, interfaces separating each domain are formed and interact only through 
coalescence dynamics: S{k) changes slowly through a self-similar process (see [17|) and the 
peak of the function will slowly move to smaller wave numbers. 



III. AN ADIABATIC ANSATZ 

Our analytic method will rely on the assumption that the intermediate region is approxi- 
mated through the growth of a periodic modulation solution of the noiseless (C-H) equation: 
we need therefore to discuss how this approach is relevant to the general case where noise is 
present and where the coarsening of the non-periodic pattern acts. Indeed, as it can be seen 
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on figure (|]), the coalescence, roughly characterized by the evolution of the position of the 
peak of the structure function, does not appear to influence the dynamics before a few hun- 
dreds of units time. At those times, the intermediate regime has ended and the modulated 
pattern is formed. More precisely, figure (0) shows the typical mean width of the pattern 
as a function of time for the same conditions as figure (Q); after a transient behavior (until 
about t = 50) where the size of the pattern is dominated by the initial conditions combined 
with the linear theory of (C-H), we observe the intermediate regime (for t between 50 and 
200 roughly). In particular, for this regime, we note that the average size of the modulation 
is Am, with a deviation of less than one percent from the value predicted by the linear theory. 
It does not mean that each modulation has a length scale of A^, but more precisely that 
the distribution of the modulations length is centered around A^, as can be seen from the 
structure factor (see figure (0)). Moreover, it suggests that the growth of each modulation 
is achieved at constant length scale, determined by the initial linear instability and thus 
centered around Am- At t around 200, the growth of the modulation is saturate (as can be 
seen from figure ([^)) and the coalescence dominates the future dynamics : the length scale 
of the structures slightly increases with time. The inset of figure (|^) shows equivalently the 
typical wave length of the modulations for the (CH) model in two spatial dimensions; it 
shows again the same plateau that is in favor of the pattern growth at a constant size. 

1.2 

1.1 
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FIG. 3. The ratio between the mean length of the modulations and the most amplified wave- 
length Am = 27r/gm as a function of time, for the same conditions than figure (^). The inset 
shows a similar curve in two space dimensions, obtained by computing the mean wave number as 
a function of time, over ten initial conditions. 

Thus, we have shown that the coalescence due to the non-periodic pattern selected at 
short time can be neglected during the growth of the modulations. We need also to quantify 
the influence of the noise during the dynamics: until now, we have simplified it to the initial 
conditions which induce then a non-periodic initial pattern which is sufficient to characterize 
the general features of the spinodal decomposition. Moreover, we have shown that the growth 
of the modulation during the intermediate regime can be considered to occur at constant 
length (centered around A^) for each modulation. But strictly speaking, noise is always 
present in the dynamics and, in addition to feeding the linear instability of the homogenous 
solution ($ = 0) for short time, it generates a systematic seed of perturbations to the quasi- 
periodic pattern. It can therefore disturb this apparent frozen dynamics at constant size. 
Figure (D characterizes its effect through the evolution of the mean length of the modulation 
for different noise levels. Each curve presents the same behavior, transient dynamics which 
selects a length scale of the order of Am, then a plateau regime (beginning around t = 50), 
which corresponds to the non linear growth, followed by a coalescence dynamics. We observe 
that the length of the plateau regime, on which we are focusing, depends strongly on the 
noise level. In particular, for low noise levels, the growth of the modulation seems to occur 
at the constant length Am, during a long period after t = 50. But, the higher the noise level, 
the shorter is this plateau. The noise stimulates the coalescence process, which interferes 
more with the pattern during the intermediate regime 0. For noise levels higher than 10~^, 
the growth regime cannot be differentiated anymore from the coalescence dynamics and in 
such cases, our assumption of growth at a constant scale would no longer be valid. 
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FIG. 4. The ratio between the mean length of the modulations and the most amplified wave- 
length Am = 27r/g function of time, for different noise level, from 10" to 10"^ (the noise 
being multiplied by 100 between each curve). The shorter the nonlinear plateau, the higher the 
noise level. Each curve is obtained through an average over 100 runs. The other characteristics of 
the simulations are the same as for figure (P). 

In conclusion, we can consider that, for low enough noise level, the non linear growth 
of the modulations is made at constant length scale and that the noise has a very weak 
influence on the dynamics of the two first stages. 

Taking advantage of this observation, we can simplify the particular study of the second 
stage of the dynamics, the non-linear saturation. The aim of this paper is therefore to 
present a detailed calculation for the growth of a periodic modulation of constant size, for 
the noiseless (C-H) equation (P in one spatial dimension for a fixed temperature Eq. 

However we hope that this approach is valid also in the limit of the small noise levels, 
where the growth of each modulation appears to be unperturbed. Although the numerical 
comparison is developed for the particular case of Am periodicity, it applies to any wave 
length. 

We can now use known results concerning non-homogeneous solutions of the Ginzburg- 
Landau equation. 
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IV. QUASI-STATIC APPROXIMATION 



In order to describe our method, we will first describe a particular family of stationary 
solutions of the one dimensional (C-H). 

For e < 0, a stationary solution exists that relies the two homogenous phases $ = ±^^^ 



\e\ 2x 

= -^tanh(^=). (3) 



Such a monotonic solution allows a continuum description of the interface between the 
two stable phases. In fact, this is a particular member of a one parameter family of stationary 
solutions of the Ginzburg-Landau equation 



-$ + 2$^ - = (4) 
These solutions, the so-called soliton- lattice solutions |T^, are : 



«Dfc,,(x) = fcASn(-, k) with e = ^ = y (5) 

where Sn(a;, k) is the Jacobian elliptic function sine-amplitude. This family of solutions 
is parametrized by e and the modulus k G [0, 1]. These solutions describe a periodic pattern 
of period 

X = 4Km, where K{k) = / ' (6) 

is the complete Jacobian elliptic integral of the first kind. This family of profiles (or alter- 
nating interfaces) can be obtain exactly as a periodic sum of single solitons and antisolitons 
(or alternating interfaces) 

^(-l)nanh(7rs(x - n)) = ^M£l^Sn(x, k) with s = and k'^ = 1 - ^ 

n K[k') 

The soliton-lattice solution can be associated with a micro phase separation locally lim- 
ited by the finite diffusion coefficient. For k = 1, Sn{x, 1) = tanh(x), we recover the usual 
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interface solution ; it is associated with a one soliton solution and corresponds to a macro- 
scopic segregation. Note that K{1) diverges ; the solution 



\e\ J\e\ 
^lA^) = -^tanh(-^x). 

is thus the limit of infinite s, when the solitons are far apart one each others (strong segre- 
gation regime). 

In the opposite limit {k — > 0, or weak segregation regime) it describes a sinusoidal mod- 
ulation 

We now seek the evolution of the solution i), according to the (C-H) noiseless 

dynamics, for a fixed reduced temperature 

{x, t) = dU^^ + 2$' - ^xx$), (7) 

Numerical simulations of that problem with a small initial condition of periodicity A, 
shows the growth of the modulation at this periodicity A. As discussed in the previous 
section, such dynamics is unstable and would, in the presence of noise for instance, loose its 
periodicity. However, we have shown that this can be neglected for low enough noise levels, 
and that this "unstable" growth of the pattern is relevant there. 

The initial condition will then be taken as the sine mode q — ^'- 

$(a;,t = 0) = p sva.{q^axx) 

where v is an arbitrary small amplitude. This profile is a member of the soliton-lattice 
family (for very small k). 

The core of the method we are using involvoes to tracking the evolution of the periodic 
modulations through a simplified equation. For that purpose, we now make the ansatz that 
at first order, these modulations belong at any time to the two parameters family of solutions 
$fe,£*, with k and e* being functions of time. Since the period is chosen to be constant and 
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equal to A, using equations and (H), we find that k and e* are related to one another 
through 

e*{k) = -2{l + e)(^^)\ (8) 

and we have eventually selected a one parameter sub-family of solutions of given spatial 
periodicity (that we will call \E'*(x, k) later on) : 

AK{k)-k^ AK{k)x 
m*{x,k) = y Sn( ^,k). 

The dynamics of ^{x,t) is now reduced to the evolution of k{t) (or equivalently e*{t)). 
Given a function $ (obtained either from experimental data or numerical simulation of 
equation (^) at time t, the ansatz assumes that there exists k so that $(x, t) ~ \l/*(x, /c). 
e*{t) can be then interpreted as a fictitious temperature : it is the temperature extracted 
from the profile at a given time, using the correspondence between e* and k of equation 
(§). For instance, at t = 0, the amplitude is small and we find that k{0) = and thus 
£:*(0) = Stt^/A^, different a priori from sq {e*{0) = y in the limit z/ — > 0, for A = Am)- In 
the same spirit, we expect that at the end of the growth, the "local temperature" of the 
interface coincides with the thermodynamic one, i.e. the quench temperature 

limt_oo£*(t) = ^0 

at which the dynmics ends. Somehow, we have assumed that the dynamics of (C-H) can 
be projected at first order onto a dynamics along the sub-family '$*{x,k), which can be 
considered as an attractor of the solutions. This is well justified when a non sinusoidal initial 
condition (of small amplitude) is chosen: we then observe, in numerical simulations, a short 
transient in the dynamics which drives the solution towards the sine mode at roughly the 
same amplitude. However, for consistency, we need to check that at any time, the solution 
of (C-H) can be well approximated by a member of the sub-family. For this purpose, we 
have developed three different algorithms, taking advantage of the general properties of 
the family of solutions ^k,e '■ either, k can be deduced both from the amplitude of the 
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oscillation equals to AkK{k)/X, or from the relation k =1- (($(A/2, t)/<l>(A/4, t))^ - 1)^ ; 
thirdly, a straightforward computation relates k to the ratio of the two first terms of the 
Fourier transform of $. We have observed that the three methods show in general similar 
results within an error of one percent. However, the validity of the ansatz has still to be 
checked by comparing the initial function t) with the extrapolated function /c) 
obtained by one of these three procedures. It is shown in figure @ at two different times 
in a numerical simulation of ; we observed that the relative differences between the two 
function is much less than 0.01. 

0.4 I ' , 1 








X 

FIG. 5. Comparison for time t = 100 and t = 140 between the numerical solution of 
(circles) and the functions ^*{x,k), with k extrapolated from the Fourier transform of <I>(x,t). 
Initial conditions are taken with v = 10~^, and A = 2-1: /qm 

Moreover, a zoom on the very early time of the one dimensional numerical simulation with 
small sinusoidal initial conditions, as presented in figure (^, shows a discrepancy between 
the different methods used to extract k from the numerics. In fact using the ratio between 
the Fourier series coefficients or the ratio between the amplitude at two specific points of 
the profiles, gives for the very beginning of the numerical simulation the value = 0, in 
agreement with the statement of the initial condition^. Indeed, these two methods account 

^Remark also that these two methods give estimations of k that are so close that they cannot be 
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for the shape of the profile, which has been taken initially as a sine mode, instead of a Sn 
mode. 

However, the use the third method, which depends on the amplitude, leads to a small 
but finite k = 0.02 due to the small but finite value of u. Even if the distortion from a 
sinusoidal function associated with this finite Jacobi modulus is small (the relative change 
in the natural period and in the shape of the function is of order 10~^), one nevertheless 
observes that the dynamics of the system is such that, within a short time, the three methods 
give results which are again in agreement. There is a short inflation period during which 
there is a change in the shape of the profile and where the shoulders of the initial sinusoidal 
inflate. That is, there exist a very short stage during which the system goes very rapidly to 
a state very close to an element of the family of the soliton lattice \l/*(a;, k). 

0.03 
0.02 
0.01 



' 2 4 6 

time 

FIG. 6. Comparison for the very early times of the dynamics between the k{t) predicted by the 
three different methods. The ratio between the Fourier series coefficients and the ratio between 
the amplitude at <I>(Am/2,t) and <I>(Am/4, t) both start at = 0. On the contrary, the third 
method, which uses the amplitude of the profile, starts from a finite value. Nevertheless, the 
three methods merge within a short time, indicating the "affinity" of the (C-H) dynamics for 
soliton-lattice solutions. 



differenciated in figure (^). 
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V. NON-LINEAR GROWTH 



Although the evokition of k{t) can be extracted from direct numerical simulations of 
as shown above, the aim of the rest of this work is to show that k{t) can be deduced via an 
explicit ordinary differential equation. Therefore, in what follows, we will seek the solution 
of ([Zp in the form 

(^{x,t) = '^*{x,k{t))+r](p{x,t) (9) 

where ip accounts for high order correction terms to \&*, while the "ansatz" assumes that 
?7 <^1 (we consider Eq and of order 1). 

To describe the evolution of the modulus k{t), or equivalently the dynamics for 6*{t), 
we will use the so-called solubility condition technique. Substituting formula (H) in the 
Cahn-Hilliard equation (|^), it gives the following dynamics 



where we have kept only the lowest order terms in the perturbation. As \E'* [x, k (t)) satisfies 
the relation : 

+ 4^*^ - 2V2^* = 0, 

we then have the following dynamics : 

dk dt 2 dx^ 'dV 'dx^^2^ ^ 

The balance of the different terms gives the small parameter of the expansion 77 ~ e* — eo ; 
we obtain ^ ^ V dt^p ~ rup. Neglecting the terms of order rj"^ in the previous equation, 
we end up solving the linear system : 

92 _ d'^* dk _ 



dx"^ dk dt dx"^ 

Here, C is the linearized C-H operator Cp = rj + 6\1'*2 — V^) p. Strictly speaking, 
this analysis is valid only for e* ; however, it is a classical assumption of the solubility 
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condition (confirmed below by the numerical results presented on figure (^) to expand it 
for the whole dynamics. 



A necessary condition for solution is that the right hand side of the system is orthogonal 
to the kernel of the adjoint operator (9^.2£)^ The Goldstone mode dk^k,e*, for £ = £* =const, 
is clearly an element of Ker(£"''), and if we consider the distribution t), such that 
^x(^)^) = ^k,e*{^), then we have d^x £Ker((9a.2£)^) |. Thus, using the scalar products 
<|> over the period A, defined as 

1 r^/s 



< f\9 >= T / f{x)g{x)d. 

A J-A/2 



X 

A/2 

we obtain the desired equation for ^ : 

n mr* dA; (£o-£*) ^ ^ *^ kK{k) (2E{k) 

(10) 



where E{k) = /g^ — k'^ sin^ x dx is the complete Jacobi elliptic integral of the second 
kind. The l.h.s can be expressed using 'ijj{x,t), defined as -^'ijj='^*, which reads : 

^(a;,t) =ln (^Dn(|, A;) - A;Cn(|, A;)^ _ lln(l - A;^). 

Cn and Dn are the Jacobi elliptic function cosine and delta amplitudes respectively. Then, 
noting that 

< dkX I dk'^k,e' >=- < dkip I dkip >= -I{k) 

where I{k) is independent of A. Finally, equation (|1^) can be recast as the following 
explicit ordinary differential equation for k{t) : 



^Remark that the partial derivative with respect to k is made with e* taken constant, since we 
are interested in a member of Ker(£^^), linearized (C-H) operator for e* . 
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In the hmit k —>■ 0, that is, for early times, equation becomes, with the wave number 
g = 27r/A associated to the period A 

Since k is proportional to the amplitude of the sine mode of wave number q, we observe 
that we retrieve the linear theory of (C-H) in that limit. 



The r.h.s. of ( ]Tl| ) is in fact proportional to 6o — £*, so that the dynamics ends when 
the fictitious temperature reaches the thermodynamic one Eq] this occurs for k = kg which 
satisfies 32(1 + kl)K{ksf = -EoX^. For A = A^, we obtain fc^ = 0.471941. This corresponds 
to the end of the non-linear growth (on figure (0 c))) and the value of kg associated with 
this steady state is well retrieved numerically by the three methods explained above. Thus, 
the asymptotic steady state solution of equation (|^) for a given period is \imt->oo^{x,t) = 
'^*{x,ks). However, no analytic solutions of equation (|lT]) have been found, and we need 
to solve it numerically. Figure (0) compares the solution of (|ll]) with the dynamic of k 
extracted by the Fourier method from the full (C-H) evolution, for the period A = A^, the 
fastest growing mode. It shows a good agreement between the two curves; in particular, 
both limit t — and t oo are well captured ; thus equation ([Tl|) remains valid even for 
k far from fc,. 
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FIG. 7. Comparison between the solution of the ordinary differential equation (|lTj) for the 
modulus k(t) (gray curve) and the modulus extracted from the full (C-H) dynamics (black line) 
with the same initial condition A;(0) = 2 • 10^"^, for A = Am- The dynamics converges to ks in both 
case for large time. The inset shows the exponential convergence of both curves in asymptotic 
regime, where the solubility condition is valid ; it compares well with kg — k{t) ~ e~^o*/^ ; in 
addition, far from /c^, the exponential growth for small time k(t) ~ e~'^o*/i6 is also retrieved by 
equation (pA]). 




VI. CONCLUSION 

We have shown that the choice of an ansatz within the sohton-lattice family allows 
a reliable description of the growth of a periodic pattern in the noiseless Cahn-Hilliard 



equation. Contrary to Ref. |19|, our ansatz relies on the hypothesis that during the first 
two stages of the dynamics, the periodicity of the order parameter remains constant. In 
this sense, it is an adiabatic ansatz. The validity of these assumptions has been investigated 
in detail and checked numerically (see Fig. |^). It enables to model the non-linear growth 
starting with spatial random initial conditions and predicts the stationary profile \E'|.^, which 
ends this non-linear growth. Although this profile might not be observable in a usual phase 



transition due to the presence of noise , we claim that this approach should be valid when 



the noise is low enough, which is the case when the quench is achieved at low temperatures. 
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We expect this approach to have a particular pertinence for axial segregation in rotating 
drums |^0[|, where the dynamics ends after the second stage. 

The use of the solubility technique combined with the choice of an adiabatic ansatz 
might be generalized to the study of other non linear dynamics. For instance, spinodal 
decomposition in superfluid Helium or Bose condensate has been argued to be described 
by a cubic-quintic non linear equation |2l]| ; in this particular case, one would first need 
to retrieve a relevant soliton-like family of solutions along which to compute the adiabatic 
dynamics. The same difficulties would arise when the method is adapted to higher space 
dimensions. 

Finally, this approach could be used to explore the self-similar scenario for coalescence 
suggested by AFM experiments for spinodal decomposition in mixtures of block co-polymers, 
depicted in Ref. [||, starting with the previous stationary distribution as initial conditions. 
The only change will be in the use of a family of solutions of growing periodicity A*, which 
would also be a slow variable of the position, since the coalescence is controlled by local 
interactions of the pattern The goal in that case would be to obtain a differential 

equation for A the order parameter with the same technique. However, these questions are 
postponed to future studies. 

The authors are grateful to David Andelman and Sergio Rica for helpful discussions. 
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